Evidence for an exotic magnetic transition in the triangular spin system FeGa2S4 
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We report positive muon spin relaxation measurements on the triangular lattice magnetic system 
FeGa2S4. A magnetic transition not previously detected by specific heat and magnetic susceptibility 
measurements is found in zero field at T* ~ 30 K. It is observed through the temperature depen- 
dencies of the signal amplitude and the spin-lattice relaxation rate. This transition is therefore not 
a conventional magnetic phase transition. Since persistent spin dynamics is observed down to 0.1 K, 
the ground state cannot be of the canonical spin-glass type, which could be suggested from hysteresis 
effects in the bulk susceptibility below Tf ~ 16 K. These results are compared to those found for 
the isomorph NiGa2S4. It is argued that the fate of the transition, which has been interpreted in 
terms of the Z2 topological transition in this latter system, is probably different in FeGa2S4. 

PACS numbers: 75.40.-s, 75.40.Gb, 76.75.-fi 
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The search for new states of matter is at the forefront 
of condensed matter research, and the geometrically frus- 
trated magnetic systems provide a fruitful playground 
[ll-Q. One of the main goals is to find and character- 
ize magnetic systems for which long-range magnetic or- 
der is absent at low temperature despite strong exchange 
interactions. Two-dimensional Heisenberg antiferromag- 
nets are good candidates since magnetic order can only 
occur at temperature T = 0. In the case of the equilat- 
eral triangular lattice, it is now believed that the system 
would exhibit 120° spin order ^4i^6:|. However, symme- 
try breaking energy terms and exchange interactions be- 
tween the layers or inside the layers beyond the nearest- 
neighbor may lead to magnetic order at finite tempera- 
ture or quantum spin-disordered ground states . Only 
a few two-dimensional triangular-lattice systems are be- 
lieved not to display magnetic long-range order: two or- 
ganic and two inorganic compounds While recent 
interesting results have been reported [9| the restricted 
amount of organic material that can be produced im- 
poses serious limitations on the measurements that can 
be performed, precluding e.g. inelastic neutron scatter- 
ing experiments. This is not the case for the inorganic 
transition metal sulfides NiGa2S4 and FeGa2S4, which 
are available in large quantities. These two insulators, 
which crystallize in the space group P3ml, consist of 
magnetic Ni^+ and Fe^"*" layers, respectively, character- 
ized by strong intralayer magnetic bonds in comparison 
to very weak inter layer couplings. The Ni and Fe ion 
sites form an equilateral triangular lattice in each layer. 

While NiGa2S4 has been studied with a wide range 
of macroscopic and microscopic experimental techniques 
fewer studies including X-ray diffraction, mag- 
netic susceptibility, resistivity, and specific heat measure- 
ments have been reported for FeGa2S4 [3, The 



available data suggest that the two magnetic systems 
are similar. They have comparable lattice parameters. 
Their Curie-Weiss temperatures 9cw are quite large, 
80 and 160 K for NiGa2S4 and FeGa2S4 respectively. 
Upon cooling below 100 K, the dc susceptibility Xdc of 
the two compounds develops weak easy-plane anisotropy 
(which is larger for FeGa2S4) but no in-plane anisotropy. 
Xdc{T) shows a bifurcation between zero- field cooled 
(ZFC) and field-cooled (FC) data at 8.5-9 K and 16 K 
(i.e. ~ for NiGa2S4 and FeGa2S4, respectively. 

The magnetic contributions to the specific heat Cm (T) of 
the two compounds are field insensitive with a depen- 
dence at low temperature, and they match one another 
after appropriate rescaling [2^. While Cm/T is negli- 
gible for T — ^ K in NiGa2S4, an appreciable value is 
measured for FeGa2S4. Finally Cm/T exhibits an un- 
usual double-peak structure: a first rounded maximum 
at 10 K for the two systems, and a second broad max- 
imum centered around 100 K and 60 K for NiGa2S4 and 
FeGa2S4, respectively. 

The microscopic magnetic properties of NiGa2S4 are 
experimentally well established. Neutron scattering ex- 
periments [lO| show that NiGa2S4 displays incommen- 
surate quasi-static short-range magnetic correlations. 
These correlations manifest themselves in a spontaneous 
field below T* ~ 9 K in muon spin rotation and relax- 
ation (/zSR) experiments [l3 . 15 1. Concerning temper- 
atures above T*: while the longitudinal-field jjSK re- 
laxation function is exponential at high temperature, it 
becomes sub-exponential or "stretched exponential" for 
T ^ T* and below. This shape indicates the presence 
of spatial inhomogencity in the relaxation rate. On the 
premise that topological Z2 vortices manifest themselves 
in the vicinity of the temperature where Cm{T) displays 
a broad maximum [2l]|, the transition at T* was tenta- 
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tively associated with the dissociation of these vortices 
In addition, the stretched exponential relaxation 
just above T* was interpreted as the signature of vor- 
tex unbinding, which is expected to give rise to magnetic 
disorder. 

We are aware of only one investigation of FeGa2S4 
by a microscopic technique. From ^''Fe Mossbauer spec- 
troscopy, a spontaneous hyperfine field has been observed 
at low temperature with an order-parameter-like temper- 
ature dependence [i^, [2^ ■ An antiferromagnetic transi- 
tion has been inferred from these data with a Neel tem- 
perature of 33 K. In this report we present /iSR mea- 
surements which are aimed at further characterising the 
nature of the transition. We also find a magnetic transi- 
tion with T* ~ 30K ~ ecw/5 rather than 9K ~ 6*01^/10 
in the nickel counterpart. Persistent spin dynamics are 
observed in both compounds down to the lowest temper- 
ature of the measurements, i.e., 0.1 K or less 2^. Never- 
theless, the two systems differ sharply in some respects. 
In contrast to NiGa2S4 no spontaneous field is detected 
in FeGa2S4, owing to a broader field distribution. In the 
latter compound no anomaly is observed in Cm{T)/T or 
Xdc(T) at T*, whereas in NiGa2S4 a weak cusp is seen 
in Xdc{T) at low fields. In addition none of the temper- 
atures of the two Cm maxima in the two systems scales 
with T* . The observation of a magnetic transition in 
FeGa2S4 while no signature for it is detected in specific 
heat and susceptibility data points to its exotic nature. 

Data were obtained from /xSR measurements car- 
ried out at the General Purpose Surface-Muon (GPS) 
(2.8 K < T < 160 K) and the Low Temperature Facihty 
(LTF) (0.10 K < T < 3.5 K) instruments of the Swiss 
Muon Source (S^iS), Paul Scherrer Institute (PSI), ViUi- 
gen, Switzerland. Most of the spectra were recorded in 
zero field, with additional spectra taken with an applied 
external longitudinal field. Measurements were mostly 
performed on a powder sample prepared in Tokyo. A 
few spectra were also recorded on a powder prepared in 
Grenoble, giving results consistent with a transition at 

/iSR techniques and their application to the study of 
magnetic materials are described elsewhere [25l - E7l |. A 
time-differential /iSR asymmetry spectrum is given by 
aor]P^^{t), where oq is the spectrometer-dependent ini- 
tial asymmetry, P^^^{t) is the muon spin polarization 
function [P°''P(t=0) = 1] [13], and ?y (0 < < 1) char- 
acterizes the so-called missing asymmetry due to rapid 
relaxation within the spectrometer dead time idt (~ 5 ns 
at PSI). This is typically caused by a broad distribu- 
tion of static fields arising from frozen magnetism. The 
remaining asymmetry arises from the muon spin compo- 
nent parallel to the static field, and its relaxation is due 
solely to dynamic processes (spin-lattice relaxation) 27 1 . 
In the following aof] will be called the effective initial 
asymmetry and aa'qP'^'^ {t) the asymmetry. 

Exploratory measurements at the GPS spectrometer 
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FIG. 1. (Color online) Examples of zero- and longitudinal- 
field /iSR spectra recorded on a powder sample of FeGa2S4. 
Solid curves: results of fits as explained in the text. The re- 
duced signal amplitude, i.e., the reduced effective initial asym- 
metry, for the 2.8 K spectrum is clearly seen graphically. Note 
that these data result from the weighted subtraction of counts 
recorded in two opposite detectors, implying that their shape 
and amplitude are not dependent of any fit parameter. 



indicated that the field history may have to be taken into 
account for T < 60 K. Hence, we first report zero-field 
measurements, for which no field was present when cool- 
ing down the sample below 60 K. In contrast to NiGa2S4, 
no early-time oscillation could be reliably resolved in the 
zero-field spectra at any temperature down to 0.1 K. 
Thus we have not detected any spontaneous mean field 
at the muon site. This and the reduced low-temperature 
value of 77, discussed below, indicate a broad distribution 
of fields. A lower bound for the width of the distribution 
at low temperature is of order AS « l/(7/itdt) — 0.2 T 
(7^ = 851.6 Mrad s^^T^^ is the muon gyromagnetic ra- 
tio). All the spectra can reasonably be described, except 
at short times [t ^ 50 ns) near T* where the static fields 
become smaller, by the stretched exponential relaxation 
function 



P'^^^it) = exp -{\zt) 



(1) 



where \z is the spin-lattice relaxation rate and /? with 
< /3 < 1 is the stretching power. Two typical zero-field 
spectra are shown in Fig. [TJ We recall that a stretched 
exponential relaxation function accounts for a continu- 
ous distribution of relaxation rates. But as shown for 
Lai_a;Caa;Mn03 this may not be justified [i^l, and a 
multi-site model might be more appropriate. Such a de- 
scription would require so many parameters (two param- 
eters per site, i.e., a weight and an exponential relaxation 
rate) that the physics would be obscured, and we shall 
not attempt it. 

Focusing first on the spectra recorded at the GPS, and 
assuming /3 to be a free parameter, we find 0.4 ^ /? ^ 
0.6 for 2.8 < T < 25 K and 0.3 < /3 < 0.4 for 25 < 
T < 35 K. The value above 40 K increases smoothly to 
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reach /? = 1 at 60 K. The final fits were carried out 
setting [3 = 1/2 somewhat arbitrarily for T < 40 K and 
leaving it free above that temperature. Indeed, when /3 
approaches 0.3, aor] becomes unrealistically large. This 
points out the inadequacy of the fitting function for t < 
50 ns. In contrast, setting /? = 1/2 yields a value for aor] 
that fairly refiects the observed signal at short times. As 
shown in Fig. [T] the resulting description of the spectra 
is reasonable. 

It is more difficult to analyze the spectra recorded at 
the LTF than at the GPS because a fraction of the muon 
beam stops in the silver cold finger of the dilution refrig- 
erator. This fraction is traditionally determined from a 
calibration measurement performed above T* [27|]. Here 
this measurement cannot be done since T* is higher than 
the maximum temperature available at the LTF. To cir- 
cumvent this difficulty spectra with overlapping temper- 
atures were recorded at both the GPS and LTF; the Xz 
values extracted from the GPS spectra were used to an- 
alyze the LTF data, with the effective initial asymmetry 
being the only free physical parameter in the fitting pro- 
cedure. With this method a reliable oq?? value was ob- 
tained from the LTF spectra. For temperatures at which 
only LTF spectra are available, this value was fixed dur- 
ing the fits. In addition, following the GPS results, /3 = 
1/2 was assumed down to 0.1 K, leaving Xz as the only 
free physical parameter. 

Figure [5] gives the temperature dependencies of aor], 
Xz and (3. We note the following three points, (i) The 
strong anomalies in aor]{T) and Xz{T) at T* ~ 30 K 
point to a magnetic transition at T*, independent of the 
details of the fitting procedure. This is a new evidence, 
after ^ ''Fe Mossbauer spectroscopy [i^l , of such a transi- 
tion in FeGa2S4. Its width is relatively large since the de- 
crease of oqt? with cooling begins at 40 K and ends only at 
« 25 K. The increase of Xz by a factor 5 on cooling from 
65 K to T* is due to the slowing down of magnetic fiuc- 
tuations often observed on cooling towards a magnetic 
transition. The drop of aor] is an indirect but perfectly 
clear indication of the onset of a spontaneous field at the 
muon site for a magnetic transition in a powder sample. 
Remarkably, there are anomalies neither in aor]{T) nor 
in Xz{T) around = 16 K where the ZFC bulk suscep- 
tibility is maximum and shows a bifurcation with the FC 
susceptibility; see Fig. [51 Conversely there is no anomaly 
in the susceptibility at T*. Furthermore, the tempera- 
ture T* does not correspond to any of the broad maxima 
observed in Cm/T. For T <t: T*, aor] is « 45% of the 
value above the transition, relatively close to the value 
1/3 that would be expected for a powder sample with 
randomly-oriented local fields. As discussed previously, 
the early-time signal relaxes rapidly and is lost in the 
spectrometer dead time, (ii) The spin-lattice relaxation 
channel is modeled by a stretched exponential function; 
the exponential relaxation expected in a homogeneous 
sample is only observed deep in the paramagnetic phase. 
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FIG. 2. (color online) The three upper panels show the 
temperature dependencies of the effective initial asymmetry 
aoT], the spin-lattice relaxation rate Xz, and the stretched- 
exponential power p measured in a powder sample of FeGa2S4 
in zero or 10 mT longitudinal fields. The solid curves for 
Xz{T) above T* result from a fit to a model explained in 
the main text. The adequacy of the model can be judged in 
the insert which displays XzjT^ versus l/T. The tempera- 
ture dependencies of the susceptibility Xdc given in CGS or SI 
units and the ratio Cm/T of the magnetic specific heat over 
the temperature are shown in the lower two panels. The 
susceptibility has been measured in an external field of 0.1 T 
after a ZFC or FC protocol. For reference dashed lines are 
plotted at 30 and 40 K. 



This suggests that the compound is magnetically inhomo- 
geneous: there is a distribution of spin-lattice relaxation 
27 , 2§- 31 1 . (iii) Since an appreciable monotonic de- 



rates 

cay of the asymmetry is observed even at extremely low 
temperature, persistent spin dynamics are present. In 
fact Xz{T) displays a plateau below « 2 K down to the 
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lowest temperature (0.1 K). This behavior seems to be a 
characteristics of frustrated mag netic compounds either 
with [32h34;] or without [s^ |36| magnetic order. Inter- 
estingly, the presence of fluctuating iron moments was 
also deduced from the Mossbauer spectroscopy study of 
FeGa2S4 



23|. 



Quantitatively, we can fit the prediction of a model 
which treats nuclear/muon spin relaxation due to criti- 
cal fluctuations of a 2D quantum frustrated Heisenberg 
antiferromagnet (37l - [39j |. to Xz{T) in the paramagnetic 
state. This model predicts for T < To/2 



Azcxr3exp(ro/r), 



(2) 



where Tq is proportional to the spin-stifFness constant 
which is related to exchange. Figure [2] shows the best fit 
to the data. We find Tq = 225 (20) and 280 (20) K for the 
zero and 10 mT longitudinal field data, respectively. The 
exchange constant J can be estimated from the relations 
keTo = Anps jsl] and Ps = [I - 0.399/(25)] J^VVS for 
5 = 2 ;37]: we find J/fcs = 8.6 (8) and 10.7 (8) K for the 
two values of To, respectively. They compare favorably 
with the value 13.3 (7) K derived from the Curie- Weiss 
temperature. 

We have also recorded spectra after field cooling the 
sample in a i?ext = 10 mT longitudinal field; see Fig. [1] 
It appears that the application of the field strongly affects 
the spectral shape at short time. This can be interpreted 
as a signature of a change of the static or quasistatic 
field distribution at the muon site. At longer times the 
spectra reflect only the muon relaxation associated to the 
dynamics of the local field. Therefore, to gain an insight 
into these dynamics we fitted the experimental data after 
truncation of the early time part, using Eq. ([1]). Cut- 
off times ranging between 100 and 300 ns were tested 
and the fitting parameters were found essentially inde- 
pendent of this cutoff. The results are shown in Fig. [2] 
The field drastically suppresses Xz, and the temperature 
dependencies of the three parameters in the vicinity of 
T* are somewhat shifted to higher temperature than in 
zero field. 

The value of i?cxt is much less than the width AB of the 
local field distribution to which the muons are submitted 
at low temperature. Therefore we would not expect the 
spectral shape at long time to be affected by Sext- This 
field effect is reminiscent of NiGa2S4 [l^ where the in- 
ternal field is much less distributed than in the current 
system, but still larger than 0.2 T at low temperature. 

Additional measurements of the field response of 
FeGa2S4 have been performed at very low temperature. 
The 10 mT spectrum measured after a ZFC protocol and 
displayed in Fig. [3] again shows an important difference 
with the zero-field spectrum. Therefore the drastic effect 
of a modest field does not depend on the temperature of 
the field application and it is still present at 0.15 K. This 
is in contrast to the field insensitivity of the specific heat 
for a field as large as 7 T [isj . 



o 



E 
E 



0.15 


1 FeGa2S4 


I I 




i 0.15 K 


• 10 mT 


0.10 


■_ ^^yw0^ 


Mi 


0.05 


iW\ - 




: 


= OmT 


0.00 




. 1 .... 1 



1 2 3 4 5 6 
Time / (lis) 

FIG. 3. (color online) Comparison of spectra measured at 
0.15 K for FeGa2S4 after zero-field cooling. The first spectrum 
is measured in zero-field and the second after the subsequent 
application of a 10 mT longitudinal field. 



The rapid increase of Xz as the sample is cooled down 
in zero field towards T* is a key result of our study. 
In a conventional magnet it would indicate a magnetic 
phase transition, see for example Ref. However, as 
no anomaly is observed in the magnetic _specific heat and 
susceptibility at that temperature [1, [l^, the detected 
anomaly in Xz does not correspond to a conventional 
phase transition. Note that a rapid increase of Xz as 
the sample is cooled down towards 9 K was also found 

[3,111]. 



for NiGa2S4 [IJ, ll5|. The transition at T* is not spin- 
glass-like, since persistent spin dynamics arc found even 
at temperatures as low as 0.1 K. 

As indicated earlier the NiGa2S4 transition at T* could 
be associated with a Z2 topological transition. In the case 
of FeGa2S4 the transition which has been revealed in the 
present study occurs at a temperature much higher than 
those of the first broad peak in the magnetic specific heat 
(«10 K) and of the bifurcation in the ZFC susceptibil- 
ity («16 K). Thus it seems unlikely that the transitions 
in FeGa2S4 and NiGa2S4 are of the same nature. This 
conclusion is further supported by the fact that hystere- 
sis effects for the latter compound are much smaller and 
occur at T*. 

In conclusion, our /iSR experiments on FeGa2S4 con- 
firm the presence of the magnetic transition at T* ~ 30 K 
in zero field which was previously observed by ^''Fe 
Mossbauer spectroscopy. With the data available we are 
unable to estimate the magnitude of the magnetic mo- 
ment involved. Because of the lack of anomalies at T* in 
the specific heat and in the susceptibility and because of a 
strong field effect, this transition is exotic. In addition, as 
usual for geometrically frustrated magnetic compounds 
and in particular for NiGa2S4, persistent spin dynamics 
are observed at low temperature. 
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